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THEORETICAL CALCULATIONS OF THE PRESSURES, FORCES, AND MOMENTS AT 
SUPERSONIC SPEEDS DUE TO VARIOUS LATERAL MOTIONS ACTING 
ON THIN ISOLATED VERTICAL TAILS 1 

By Kenneth Margolis and Percy J. Bobbitt 


SUMMARY 

Velocity potentials, pressure distributions, and stability de- 
rivatives are derived by use of supersonic linearized theory for 
families of thin isolated vertical tails performing steady rolling, 
steady yawing, and constant-lateral-acceleration motions. Ver- 
tical-tail families ( half-delta and rectangular plan forms) are 
considered for a broad Mach number range. Also considered 
are the vertical tail with arbitrary sweepback and taper ratio at 
Mach numbers for which both the leading edge and trailing edge 
of the tail are supersonic and the triangular vertical tail with a 
subsonic leading edge and a supersonic trailing edge. For pur- 
poses of completeness, analogous expressions and derivatives for 
sideslip motion obtained primarily from other sources are 
included. 

Expressions for potentials, pressures, and stability derivatives 
are tabulated. Curves which determine the stability derivatives 
for half-delta and rectangular tails are presented which enable 
rapid estimation of their values for given values of aspect ratio 
and Mach number. In order to indicate the importance of end- 
plate effects, several comparisons are shown of the derived results 
(based on a zero-end-plate analysis ) with those corresponding to 
a complete-end-plate analysis. 

INTRODUCTION 

Detailed knowledge of the loading, forces, and moments 
acting on vertical tails undergoing various maneuvers is a 
necessary prerequisite for determining the lateral dynamic 
behavior of aircraft. The information presently available is, 
in many instances, insufficient to enable reliable estimates to 
bo made of the vertical-tail contributions to airplane stabil- 
ity at supersonic speeds. Aside from calculations for several 
“slender” configurations, theoretical results to date have 
been concerned, for the most part, with tail configurations 
either subject to a constant sideslip attitude or performing 
a steady rolling motion (see refs. 1 to 13). 

For the sideslip motion, the effects of Mach number and 
aspect ratio on the aerodynamic loads of a number of tail 
configurations with both one and two planes of cross-sec- 
tional symmetry have already been investigated extensively. 
The same effects on tail arrangements in a rolling motion 


have also received considerable attention which has mainly 
been directed toward tails with two planes of symmetry such 
as cruciform arrangements. Additional theoretical analysis 
devoted to the evaluation of the Mach number and aspect- 
ratio effects on the forces and moments acting on tail sys- 
tems in roll with one plane of cross-sectional symmetry is 
required. Tail arrangements performing a steady yawing 
motion or a constant-lateral-acceleration motion have re- 
ceived little attention to date in the literature. Yet the 
forces and moments produced by these motions are by no 
means negligible, and some indication of their magnitudes 
is necessary, particularly at supersonic speeds, in order to 
evaluate their relative importance on lateral stability. 

The primary purpose of this report is to present the results 
of a theoretical investigation to determine at supersonic 
speeds the pressures, forces, and moments acting on several 
families of thin isolated vertical tails subject to various 
lateral disturbances. Three motions axe treated: steady roll- 
ing, steady yawing, and constant lateral acceleration. A 
fourth motion, namely, constant sideslip, although analyzed 
previously, is included for purposes of completeness. The 
basic plan forms considered are: (a) half-delta tail with 
either a subsonic or supersonic loading edge, (b) rectangular 
tail, (c) general sweptback tail of arbitrary taper ratio with 
supersonic leading and trailing edges, and (d) triangular tail 
with a subsonic leading edge and a supersonic trailing edge. 
The half-delta and rectangular vertical tails are analyzed in 
detail in that forces and moments (expressed in the form of 
stability derivatives) and their variations with Mach num- 
ber and aspect ratio are presented in a series of simple charts. 
Useful expressions and formulas are included for the other 
plan forms which enable similar calculations to be carried out. 

A secondary objective, in view of the geometric nonplanar 
characteristics of tail arrangements, is to consider the esti- 
mation of the mutual aerodynamic interference that exists 
between the vertical and horizontal tails. In order to gain 
some insight into tbe possible effects of such interference, 
several of the derived results are compared with correspond- 
ing calculations based on a complete-end-plate analysis. (A 
complete-end-plate analysis implies that the horizontal tail 
acts as a perfect reflection plane.) 


’Supersedes NAOA Technical Notes 3373 by Kenneth Margolis, 1966, and 3S40 by Peroy J. Bobbitt, 1954. 
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SYMBOLS 
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<r=m 
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B=-jW=i 

p,r 

0 
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p 

2 

AP/g 


<P 


A<p 


<Pd 

<Pl 

A(x,z) 

Cr 

b 

s 

X 


A 

m 


rectangular coordinates used in analysis (see 
fig. 2(a)) 


rectangular coordinates of doublet 

longitudinal and vertical distances, respectively, 
that origin is displaced relative to tail apex 
(see fig. 2(c)) 

perturbation velocities in x- and ^/-directions, 
respectively 

free-stream or flight velocity (see fig. 2) 

Mach number, T/Speed of sound 


angular velocities about x- and 2 -axes, respec- 
tively (see fig. 2) 
angle of sideslip (see table I) 
rate of change of 0 with time, dfi/dt 
time 

density of air 


free-stream dynamic pressure, ~ pV 2 

£ 


coefficient of pressure difference between op- 
posite sides of tail surface due to particular 
motion under consideration, positive in sense 
of positive side force (see fig. 2) 
constant determining degree of homogeneity of 
quasi-conical velocity field 
perturbation velocity potential due to particu- 
lar motion under consideration evaluated on 
positive y-side of tail surface (see fig. 2) 
difference in perturbation velocity potential 
between two sides of tail surface, 
( p( x fi + , z )—<p(X)0~,z) 

potential of supersonic doublet distribution 

potential of line of doublets 

doublet-strength function 

line-doublet-distribution function 

root chord of vertical tail 

span of vertical tail 

area of vertical tail 

taper ratio of vertical tail, Tip chord/Root 
chord 

aspect ratio of vertical tail, b 2 /S 
slope of leading edge of tail; cotangent of 
sweopback angle of leading edge (see fig.l) 


1— A /1 
Bin 


E(Bm) 


■J^L-^ll-BBn 2 ) 

E’{k) 


E’{k) 


complete elliptic integral of second kind with 
modulus ■jl—k 2 , 



(1— k*) sin 2 n dn 


K'(k ) complete elliptic integral of first kind with 

modulus Vl— k 2 , 

C T/2 dn 

Jo yi — (1 — k*) sin 2 n 

r, cc, to arbitrary constants (co=<aBm) 

, ^fl+F[k 2 (.l+k r )K / (k) + (l-4:k 2 +k i )E'(k)] 

Tp (2 — F) (l—2k*)E'(ky+k 3 (l+ ¥)K' (k) E' (Ic) - k<K' (k ) 2 

/= (l+E) w [(l+k 2 )E'(k)-2lc‘K'(k) ] 

2[(2-E)(l-2E)E'(kf+k 2 (l +k*)K' (k) E' (k) —k*K' (kf\ 


k^Jl+lPiil+l^E'W^lPK'ik)} 

Tr (2 -IP) (I-2k*)E'(ky+k i (l + t*)K'(k)E' (Jc) -k*K'(ky 

, (l+^Vil-lP+k^E'W-Pil+k^K'm 

“ r 2k{(2-JP)(l-2k*)E'(ky+m +Jc r )K' (k) E' (1c) — k*K' (A*) 2 ] 

e infinitesimally small quantity 

x—B?z<t 

7 ~ Vl-BV yP-REif+A) 

Yo expression for 7 indicated at y= 0, — == - L . 

y ’ yi-svyi -bw 

4>,v variables used for integrating purposes 

Y side force 

N yawing moment 

L' rolling moment 

C y side-force coefficient, Y/gS 

C H y awing-momen t coefficient, N/gSb 

Ci rolling-moment coefficient, L'/gSb 
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Subscripts: 

p refers to rolling condition 

r refers to yawing condition 

r=l refers to unit-yawing condition 

/3 refers to sideslip condition 

/3=1 refers to unit-sideslip condition 

/S refers to constant-lateral-acceleration condition 

w based on wing dimensions 

1,2 components used for fi derivatives 

Abbreviations : 

L. E. leading edge 
T. E. trailing edge 

All angles are measured in radians. 

ANALYSIS 

SCOPE 

The vertical-tail plan forms considered herein are shown 
in figure 1. The leading edge has arbitrary sweepback and 
the trailing edge may be either sweptback or sweptforward. 
The permissible ranges of sweep, aspect ratio, and taper 
ratio for the supersonic-leading-edge configurations are 
determined by the conditions (indicated in fig. 1) that both 
the leading edges and trailing edges remain supersonic and 
that the Mach line emanating from the leading edge of the 
root chord does not intersect the tip chord. Also, the tip 
chord must be parallel to the root chord. For the subsonic 
leading-edge configurations, only the case of zero taper ratio 
is considered and the restriction to supersonic trailing edge 
is imposed. 

Expressions based on linearized supersonic-flow theory 
are obtained for the perturbation velocity potentials and 
pressure distributions due to steady rolling, steady yawing, 
and constant lateral acceleration. For purposes of complete- 
ness, analogous results for constant sideslip motion obtain- 
able, in general, from references 9 and 10 are included. 
The expressions, which are derived for the condition of zero 
geometric angle of attack and which are valid for low rates 
of angular velocities, small sidelip angles, and small angle- 
of-sideslip variation with time, are tabulated so that they 
may be utilized conveniently in the calculation of load 
distributions and the corresponding forces and moments. 

Two important members of the family of vertical-tail 
plan forms are considered in detail. These are the rectangu- 
lar tail and the triangular tail with an unswept trailing edge, 
that is, the half-delta tail. For these tails, closed-form 
expressions are derived for the side force, yawing moment, 
and rolling moment due to each motion. The resulting 
formulas are expressed in the form of stability derivatives 
and are tabulated; simple charts are presented which permit 
rapid estimation of the 12 stability derivatives for given 
values of aspect ratio and Mach number. Tabulation 


of the derivatives for subsonic-edge triangular tails with 
trailing-edge sweep are also presented. 

Three systems of body axes are employed in the present 
report. For plan-form integrations and in the derivation 
and presentation of velocity potentials and pressures, the 
conventional analysis system shown in figure 2 (a) is utilized. 
In order to maintain the usual stability system of positive 
forces and moments, the axes systems shown in figures 2 
(b) and 2 (c) are used in formulating the stability derivatives. 
A table of transformation formulas is provided which enables 
the stability derivatives, presented herein with reference to 
a center of gravity (origin) located at the leading edge of the 
root chord (fig. 2(b)), to be obtained with reference to an 
arbitrary center-of-gravity location (fig. 2 (c)). 

BASIC CONSIDEBATION8 

The calculation of forces acting on the vertical tail essen- 
tially requires a knowledge of the distribution of the pressure 
difference between the two sides of the tail surface. This 
pressure-difference distribution is expressible in terms of the 
perturbation-velocity-potential difference or “potential jump 
across the surface” A<p by means of the linearized relationship 

( V ^ A<p+ ^i A<p ) (1) 

Inasmuch as for the present investigation thin isolated tail 
surfaces are considered and thus no induced effects are 
present from any neighboring surface, the perturbation 
velocity potentials on the two sides of the tail are equal in 
magnitude but are of opposite sign. Equation (1) may then 
berewritten in terms of theperturbation velocity potential <p as 
follows: 

*£=i_ ( M 

where <p is evaluated on the positive y-side of the tail surface. 

The basic problem, then, is to find for each motion under 
consideration the perturbation-velocity-potential function 
<p for the various tail regional areas formed either by plan- 
form or plan-form and Mach line boundaries. (See, for 
example, the sketch given in table I.) 

For time-independent motions, such as steady rolling, 
steady yawing, and constant sideslip, the potential functions 
are ' of course independent of time (i. e., the last term in 
eqs. (1) and (2) vanishes) and may be determined for the 
subsonic-leading-edge cases by the doublet-distribution 
method of references 14 and 15. The details of the method 
and its application are given in the appendixes. The super- 
sonic-leading edge configurations are analyzed by the well- 
known source-distribution method utilizing the area-can- 
cellation — Mach line reflection technique of reference 16. 
The mathematical details are not presented herein, because 
it is felt that previous papers dealing with wing problems 
(e. g., refs. 17 to 20) have applied the basic method in suffi- 
cient detail. The main difference to be noted is that the 
root chord of the isolated vertical tail is, in effect, another 
free subsonic edge similar to the tip chord and must be 
treated accordingly. Actually, tail regions I and III (refer 
to the sketch in table I) are not affected by the additional 
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tip, and wing results in these regions for constant angle 
of attack (ref. 20), steady rolling (ref. 20), and steady pitch- 
ing (ref. 18) are applicable to constant sideslip, steady 
rolling, and steady yawing motions, respectively, for the 
vertical tail, provided appropriate changes in coordinates 
are introduced and the proper sign convention is maintained. 

The time-dependent motion considered in the present 
report, that is, constant lateral acceleration, can be analyzed 
in a manner analogous to that used for a wing surface under- 
going constant vertical acceleration (e. g., refs. 21 to 23). 
By following this procedure, the basic expressions for the 
perturbation velocity potential and pressure coefficient 
(evaluated at time t= 0) may be derived as follows: 


<p= 

- 

1 

U 1 

(3) 

A P P f 

~ M * ( — ^ ■ 

M*x /AF\ ,4 


2 " 5*1 

V 2 )r - 1 

f V \ 2 

W 


Equations (3) and (4) may be deduced directly from equa- 
tions (1) and (2) of reference 23, provided the corresponding 
tail motion is substituted for each wing motion and, further, 
that care is exercised in preserving the conventional system 
of positive forces and moments. The choice of time t=0 in 
equation (4) was made for pin-poses of convenience and 
simplicity, inasmuch as the pressure due to constant sideslip 
is eliminated, and thus only the increment in pressure due 
to time rate of change of sideslip, that is, p, remains. 

The right-hand sides of equations (3) and (4) are composed 
of terms involving steady or time-independent motions, in 
particular, the motions previously discussed in this section. 
Thus, once the potentials and pressures are determined for 
steady yawing and constant sideslip, corresponding ex- 
pressions may be obtained for constant lateral acceleration 
by use of equations (3) and (4). 

Derivations of the potentials and pressures for the various 
regions of the tail plan forms under consideration have been 
carried out for each motion by using the methods and tech- 
niques discussed. Tabulations of the potential and pressure- 
distribution functions are given in tables I and EE for con- 
stant sideslip, in tables m and IV for steady rolling, in 
tables V and VI for steady yawing, and in tables VH and 
VLLL for constant lateral acceleration. 

The forces and moments acting on the vertical tail due to 
each motion may be obtained by plan-form integrations of 
the appropriate potential and pressure functions and may 
be given as follows (the center of gravity is assumed to be at 
the leading edge of the root section) : 


fTiP f T. E . A P 

q — ax az 

J Root Jle. 2 

(5) 

fTlP fTE-AP , , 

-q — x ax az 

J Root JL.E. 2 

(6) 

r tip r t. e. ap 
q 1 — z ax az 

J Ro.t J L. E. 2 

(7) 


For steady motions, ^-=^> and thus the first integration 

with respect to a; in equations (5) and (7) yields <p; hence, 
equations (5) and (7), when applied to steady motions, reduce 


to essentially a single integration involving the potential 
function. 

The nondimensional force and moment coefficients* and 
corresponding stability derivatives are directly obtainable 
from the definitions given in the list of symbols. For 
example, 



Inasmuch as the various pressure coefficients are linear with 
reference to their respective angular velocities, attitude, or 
acceleration (i. e., linear in p, r, p, or p), the partial deriva- 


tive in the preceding example may be replaced by 
and the derivative C'„ r is then expressed as 


1 

rb/V 


1 fTiP r t.e. A p 

C„ r = —xdxdz 

ail r jBootjL.E. 2 


sb v 


( 8 ) 


Corresponding expressions for the 11 other derivatives 
C«p C lfl , C Yr , C h , C Yp , C nf , Ci p , C r . , (7 B , and C,^ may bo 
obtained in an analogous manner. 

In the present report, the triangular vortical tail with 
unswept trailing edge (half-delta) and the rectangular vor- 
tical tail have been analyzed in detail. The results obtained 
upon performing the plan-form integrations and other 
mathematical operations indicated in the previous discussion 
are tabulated in table IX. Table X presents similar results 
for the subsonic-edge triangular tails with trailing-edge sweep. 
For convenience, a table of transformation formulas is pre- 
sented (table XI) which enables the derived results for the 
stability derivatives (tables IX and X) to be expressed with 
respect to an arbitrary center-of-gravity location. 

Values of the elliptic-function parameters appearing in the 
subsonic-edge formulas are presented graphically in figure 3. 


COMPUTATIONAL RESULTS AND DISCUSSION 
The formulas for the stability derivatives given in tablos 
IX and X are seen to be functions of t he tail-a spect ratio A 
and the Mach number parameter B=^JM 2 — 1. Use of the 
combined parameter AB for the abscissa variable and appro- 
priate choice of derivative parameters for the ordinates allow 
the analytical results for most of the stability derivatives to 
be expressed graphically by means of a single simple curve; 
the stability derivatives due to /3-motion require two curves. 
Figures 4 to 9 present the results for the half-delta tail and 
figures 10 to 14 present the analogous results for the rec- 
tangular tail. The lower limit AB= 1 for the rectangular 
vertical tail corresponds to the condition where the Mach 
line from the leading edge of the root chord intersects the 
trailing edge of the tip chord. Values of the derivatives for 
the situation where the Mach line from the leading edge of 
the root chord intersects the tip chord, that is, values of 
AZ?<1, cannot, in general, be obtained easily because of the 



THEORETICAL AERODYNAMICS OF THIN ISOLATED VERTICAL T AILS AT SUPERSONIC SPEEDS 


389 


fact that the calculation involves the consideration of inter- 
acting external flow fields. The lower limit AJB=1 in this 
case is not very restrictive except for very low aspect ratios 
at low supersonic Mach numbers. 

In considering the curves given in figures 4 to 14, the fol- 
lowing facts should be kept in mind: (a) The results are for 
a completely isolated vertical tail, (b) the center-of-gravity 
location is assumed to be at the tail apex, and (c) parameters 
used for nondimensionalizing purposes are the area and span 
of the vertical tail. For other center-of-gravity locations, 
analogous curves may be drawn by use of the presented data 
and the axes-transformation formulas given in table XI. 

Thus far, only the isolated Vertical tail has been consid- 
ered, that is, the zero-end-plate solution. For comparison 
purposes, results for several of the derivatives based on a 
complete-end-plate analysis have been obtained and are 
presented in figures 15 and 16. The comparisons are shown 
for the side-force and yawing-moment derivatives due to 
constant sideslip, steady yawing, and constant lateral accel- 
eration. The complete-end-plate results for these tail de- 
rivatives are obtainable from stability derivatives previously 
reported for symmetrical wings (refs. 18, 20, 23, and 24), 
provided modifications are introduced to account for (a) 
changes in nondimensionalizing parameters, (b) correspond- 
ence of wing and vertical-tail motions, and (c) preservation 
of sign convention for positive sense of motion, moments, 
and so forth. The transformations of symmetrical-wing 
derivatives into complete-end-plate derivatives for vertical 
tails having the same plan-form geometry as the half-wing 
may be summarized as follows : 

C Yfi = — 1 X (Expression for C La with wing aspect ratio re- 
placed by 2/1) 

^ 4 1+A+X\ 




C: 


3 A (1+X) 2 
ratio replaced by 2/1) 
2 1+X+X 2 


X (Expression for C m with 'wing aspect 


X (Expression for C L with wing aspect 


3/1 (1 + X) 2 
ratio replaced by 2 A) 

Cn=^ 2 - 1 q^) X 4 ~ X (Expression for C Mq with wing aspect 

ratio replaced by 2 A) 

2 1+X+X 2 

(1+X) 2 X (Expression for C h - a with wing aspect 


ratio replaced by 2A) 

C ^ = ~9Z 2 ° d+tr ) X (Expression for G nk with wing 
aspect ratio replaced by 2/1) 

where C La , C„ a , C Lq , C„ q , C L . a , and C n - a are conventionally 
defined wing derivatives (see, for example, ref. 23). Figures 
15 and 16 indicate that for a given aspect ratio the effect of 
an end plate decreases with increasing Mach number and 
that for a given Mach number the effect of an end plate 
decreases with increasing aspect ratio. Although these con- 
clusions apply specifically to those derivatives presented, it 
is felt that similar evidence would be found for the other 
derivatives as well. The percentage differences between 
zero- and complete-end-plate results vary of course with 
center-of-gravity location, as well as with Mach number and 


aspect ratio, but in general are not too large for the side-force 
and yawing-moment derivatives considered except at the 
lower values of AJB. 

The stability derivatives as presented herein have been 
made nondimensional with respect to vertical-tail parameters 
such as tail span b, tail area S, and the angles pb/V, rb/V, 
and fibjV. The magnitudes of the derivatives may, there- 
fore, appear to be quite large with respect to the expected 
tail contributions to the derivatives for a complete airplane. 
The following factors should be used in converting the pre- 
sented analytical and numerical results to corresponding 
derivatives (denoted in the following relationships by sub- 
script w) based on wing area S a , wing span b w , and angles 
pb„l2V, rbj2V, and fibj2V: 

( Cr *)rK ° 7 * 

(?*)*’ K ° ! e) 

( 0 y v)j ( Cr ): (° Y b)r 2 £ h ( Cy *' Cy " Cr b) 

CONCLUDING REMARKS 

Velocity potentials, pressure distributions, and stability 
derivatives have been derived by use of supersonic linearized 
theory for families of thin isolated vertical tails performing 
steady rolling, steady yawing, and constant-lateral-accelera- 
tion motions. Vertical-tail families (half-delta and rectangular 
plan forms) are considered for a broad Mach number range. 
Also considered are the vertical tail with arbitrary sweepback 
and taper ratio at Mach numbers for which both the leading 
edge and trailing edge of the tail are supersonic and the 
triangular vertical tail with a subsonic leading edge and a 
supersonic trailing edge. For purposes of completeness, 
analogous expressions and derivatives for sideslip motion 
obtained primarily from other sources are included. 

The effects of a complete end plate on several of the side- 
force and yawing-moment derivatives have been considered, 
and it appears that only for relatively smal l value s of the 
aspect-ratio — Mach number parameter A^JM 2 —1 do the 
complete-end-plate and zero-end-plate values differ signifi- 
cantly enough to warrant further study of finite-end-plate 
corrections. 

An additional point of interest pertinent to the present 

investigation is that the results obtained for the yawing- 

moment derivatives due to steady yawing and constant 

lateral acceleration <7. and O n may be used to approximate 

6 

the aerodynamic damping of the lateral oscillation in yaw to 
the first order in frequency. This approximation to the 
lateral damping is given by the expression C„ r — C n and can 

P 

be rapidly calculated from the curves and formulas given 
herein. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., March 5, 1956. 
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APPENDIX A 


DETERMINATION OF PRESSURE-DISTRIBUTION EXPRESSIONS FOR A SUBSONIC-EDGE TAIL UNDERGOING YAWING AND 

ROLLING MOTIONS 


A method for solving supersonic-flow boundary-value 
problems governed by the classical, linearized, partial- 
differential equation 


jp &p &<p 

dx s by 3 bz 3 


(Al) 


has been developed in reference 14 and an application to 
rolling and pitching triangular wings is given in reference 15. 
This method allows the prediction of the disturbance- 
potential function <p, and hence the pressure distribution, for 
planar lifting surfaces. The analysis given in reference 15 is 
briefly summarized herein and is applied to the determination 
of the pressure distributions on a triangular-vertical-tail 
surface (fig. 2) performing rolling and yawing motions. 
(Tawing in the xz-plane is analogous to pitching in the ay- 
plane.) 


DETERMINATION OP THE FORM OF THE VELOCITY POTENTIAL 

As is well-known, the potentials of both the supersonic 
source and the supersonic doublet and their distributions 
represent solutions of equation (Al). For the determination 
of the potentials and pressure distributions of lifting surfaces 
with subsonic leading edges, a distribution of doublets that 
uniquely satisfies the prescribed boundary conditions is 
required. The boundary conditions on the vertical tail 
for the motions to be considered herein are as follows: 

On the rolling vertical tail, 


2 

v=pz=xp-=xpQ (A2) 

and on the yawing vertical tail, 

v——rx (A3) 

In addition, the following relations must-be valid on the 
surfaces of the tail: 

For the rolling motion, 



and for the yawing motion, 



It is also necessary that the pressure along the stroamwiso 
edge be zero. 

The potential in space produced by a distribution of 
doublets, for example, in the xz-piane, with the doublet 
axes normal to the plane is 


* 2) -| JI vg =gg w§= 0 =m m 


where the area S is the region of the xz-plane intercepted 
by the forecone from the field point (x, y, z). 

The potential on the surface carrying the doublet distribu- 
tion is given by 


tpD (x,z) — Ivm 

»->±o 


["A rr -Afe m ^ i 

jJs v(x-?) a -5 i (z-r) j -5v_ 


As stated in reference 15, this surface potential is directly 
proportional to the doublet-strength function A(x,z); that is, 


<Pd(x,z),.± 0 = ± vA{x,z) 

The surface-pressure velocity u(x,z) then becomes 


and the linearized pressure coefficient 

A P 4w(x,z) y _ + o 
2 ~~ V 

may be written as 

A P 4ir bA(x,z) 

y V 5x 


(AO) 


(A10) 


(All) 
(A 12) 


The problem to be considered in this appendix is one in 
which the sidewash on tho surface is prescribed (see oqs. 
(A2) and (A3)) and the surface velocity potential has to be 
determined. The doublet-strength function A(x,z) then is 
an unknown and the determination of this quantity requires 
in general the solution of an integral equation. In some 
cases the general form of the surface-potential function 
A(x,z) is known or can be obtained by inverting an integral 
equation. The problem then resolves simply into an evalu- 
ation of the arbitrary constants of the general solution by 
making use of the prescribed boundary conditions. 

Brown and Adams in their analysis of rolling and pitching 
triangular wings with subsonic leading edges (ref. 16) were 
able to determine the function A(x,z) by utilizing the concept 
that the conical properties of the produced flow gave rise to 
potentials and pressures in the crossflow planes that wore 
similar in form to the potentials and pressures acting on flat 
finite segments in a two-dimensional flow; these segments 
correspond to a section of the wing in any crossflow plane. 
This remarkable connection between linearized supersonic 
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conical flow and incompressible two-dimensional flow is 
discussed by Busemann in reference 25. 

A more general and rigorous approach to obtain the 
doublet-strength function may be formulated from an 
analysis presented in a later paper by Lomax and Heaslet 
(ref. 26) dealing also with conical and the so-called quasi- 
conical problems. In this analysis a general surface-pressure- 
coeflicient expression has been dete rmin ed for planar lift ing 
surfaces with prescribed boundary conditions of the form 

v~z‘ 9 (£) (A 13) 

This expression is 


The arbitrary constants r and to in the so-called distribution 
function /(0) are, in terms of bi and b 3 , 


V ( 2&i . b 3 \ 
27 tB \3m 2 ‘ m / 

V 6, 

: 2t tB 3 m 2 


By relating equation (A16) to equation (A9), the doublet- 
strength function A (x,z) is seen to be 



(A18) 


AP = /*V 

2 {By t-o V(m— d)(d— mi) 


(Al4) 


where b f are constants, 0 =— > k is determined by the boundary- 

condition equation (eq. (A13)), and m and are the tan- 
gents of the apex angles of the two panels of the lifting sur- 
face. When rrii=m, the lifting surface is symmetrical about 
the common root chord of the two panels, and when m^m, 
the lifting surface is asymmetrical about this chord. From 
equations (A12), (A14), and (A9) the form of A(x,z) or, 
synonymously, the form of the surface potential may be 
obtained by a simple integration. It should be mentioned 
at this point that reference 26 presents a method for deriving 
the arbitrary constants b t in the pressure coefficient (eq. 
(A14)). This method is related to that of reference 15 
which concerns itself with obtaining the arbitrary constants 
in the velocity potential. 

By application of equation (A14) to the boundary problem 
of the triangular vertical tail (mi=0) and by noting from the 
prescribed boundary conditions (see eqs. (A2) and (A3)) that 
k= 1, the pressure coefficient for both the yawing and rolling 
motions is 


A P x &o+&i0+M 2 /i-rs 

T ~b~W=W ’ 

The constant b 0 in this expression must be set equal to zero 
in order to satisfy the condition that along the streamwise 
edge the pressure must be zero. 

The velocity potential on the vertical-tail surface is easily 
obtainable from the pressure expression by the formula 


v r ap ^ 

V~~t - — &£ 

4 Jl.b. 2 


and has been found to be 


<p = TTX ! f(- 


where 


j (~)=/( 0 ) = i-rO+um) V0(m— 0) 


(A 16) 
(A 17) 


A comparison of the potential given by equations (A16) 
and (Al7) and the potential obtained for the slender, rolling, 
vertical tail reported in reference 5 shows, as expected, that 
both are of the same form. 


EVALUATION OF THE CONSTANTS r AND o, 


The constants r and &> in the expressions for the velocity 
potential given by equations (A16) and (A17) are still to be 
determined. As indicated previously, the expression for 
the pressure coefficient, and hence the velocity potential, can 
be determined completely through an application of the pro- 
cedures developed in reference 26; however, many of the 
integrations and integrating procedures required in the 
method in reference 15 were already known to the junior 
author at the inception of this project and, for this reason, 
the analysis herein to determine the constants r and a> 
closely parallels the procedures discussed in reference 15. 

The determination of the constants r and &> depends upon 
satisfying the boundary conditions associated with the verti- 
cal tail for the rolling and yawing motions. These boundary 
conditions are given by equations (A2) to (A7). The needed 
expressions for the prescribed velocities and their derivatives 
with respect to 0 in terms of the distribution function /(</) 
are derived in appendix B. 

For the rolling motion the constants r and o> may be 
obtained by replacing f(a) by its equivalent (eq. (A17)) in 

b(v/x) 

the equations given in appendix B for v/x and and then 


applying the boundary conditions given by equations (A2) 
and (A4). When the integrations have been performed, 
the resulting equations may be solved simultaneously for 
t„ and cop. The yawing constants are obtained in a like 
manner with equations (A3) and (A6) replacing equations 
(A2) and (A4), respectively. 


In the calculation of the quantities vjx and 


<Kg/g) 

dd 


> any 


value of 8 may be considered. It is advantageous for inte- 
gration purposes to let this value of 8 be zero. However, 
since one of the limits of integration is zero and since in 
several of the integrands a singular point exists at 0= cr=0, the 
integrations in which these singularities occur must be per- 
formed for 0 arbitrary before 0 can be set equal to zero. 
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Substituting equation (A17) into equation (B4) gives, for 9 equal to zero, 


» = _i ~ 

x Bj o _ 


— 3(r5tr+&jBm) — Bo) tanh 1 Vl — B 2 a 2 . 2(rB<r+aj5m) -jBaiBm — Ba) 

(1 -BV) 5/2 H MW . 


d(5«r) + 


sfe BT* 

2 (tBB+coBiti) -jB9(Bm—B9) t/T—B*(P 


1 f fim (rgff+ggwQ s/BajBm-Bo-) (1--B+0) 2 Vl-W 




(1 -B?<?y{B<x-Bey diB<x) ' 


B*6 


} 


Carrying out the integrations in equation (A19) yields 


* B-Jl+Fil-k 3 ) 1 - 


{S , (4)[24 3 T+£4 2 (l+4 i )]-S'(4)[T4(l+4 s )-i3(l-44 2 +4 4 )]} 


where 


and 


k= 


1— Vl —BW 


Bm 


(A10) 

(A20) 

(A21) 


j=ioBm=c 


24 

1+4* 


These integrations were accomplished with the aid of the tables in references 27 and 28 and are discussed in appendix C. 
Substituting the distribution function into equation (B5) results in, for 9 approaching zero, 


d(g/s ) 

b9 


rBmr. 

Jo 


lim 

e-*o 


‘iBairBa^r aBrri) s]Ba{Tlm — Ba) tanh ‘Vl — BV 3 Ba(rBa-\-ojBm) -<jBa(Bm — Ba) 
(1 -B'o*)*'* (1 -£V) 2 . 

Pl30(r7?g+ai7?m) ^jBa(JBm — Bo) ( rBa-ycoBm ) -jBa(Bm — Ba) , 


d{Ba) + 


feU 


L Vi -w QBa-Bey 


Vl -w (l-B*o*)(Bv-Be) 


2(jBa+ utBrri) , JBa{Bm-Ba ) Vl-W 
(Ba-BOf 


d(Ba) + 


r p 

Jb(9+.)L 


B9(rBa-\-CiiBm) -yJB a (Bm, — jBtr) 


(r5<r+col?m) ■yjBa(Bm — Ba) , 2(rBa-{-oiBm) ■jBa^Bm—Ba) V 1 — B^O 3 
Vl -2W (1 (Bcr-BO) + {Ba-Bey . 


Vi-w (Ba-Bey 
d(Ba) — 


2B0(tB0-{- oiBm) -JB9{B9^ iBtri) _2[— iB i 9 2 T-\-B9Bm (3r — 2a>) + uBhn*] -JT—Wd * ^ 
JJeVl-W BejB9(B9-Bm) J / 

By performing the integrations in equation (A22), the following expression is obtained: 


b(vjx)_ 


[k*K'(k) (t+t4 2 +2S4) + S'(4) (2t4*-2t-54-2t4 4 -54 3 )] 


vt+F(i-4*)* 

Consider the rolling case for which t=t„ and co=co p . Also from equations (A2) and (A4), for 0=0, v/x=0 and 

- 24 

Solving equations (A20) and (A23)for t„ and co p , with co=B7na) P =Yq-^jW p , gives 


(A22) 

(A23) 
b(v/x) 


50 


—p -JT+F [4*(1 ■ +4*)X' (4) ■ + (1 — 44*+4 4 ) S' (4)] 

Tp 7r[-4 4 S'(4) 2 +4 2 (l+4 2 )S , (4)S'(4) + (2-P)(l-24 2 )S , (4) 2 ] 

p(l+4 s ) 3 *[24 I K'(4)-(l+4 s )S'(4)] 
“ p_ 27r[-4 4 s:' (4) J +4* (1 +4 i )S' (4) S' (4) + (2— 4 s ) (1 -2k 2 ) E' (4)^ 


-=p. 

(A24) 

(A26) 
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For the yawing case r=T r , u r , and, from equations (A3) and (A6) ^=— r and 

equations (A20) and (A23) simultaneously after making these substitutions yields 


=0, respectively. Solving 


Brk^l+k l [2lc?K , {k)-(l+lc ! )E'{k)] 

Tr {ky+V(l+W)K' Qc)E' (jfe) + (2-F) (1-2 ¥)E' (fc) 2 ] 

-Br (1- +P) 3/1 [P(1 +V)K' Qc) -2 (1 -&+¥) E' (Jc)] 

“ r 2Trk[-k i K'(ky+k a (l+k*)K , (k)E'(k)+(2-k 2 )(l-2k*)E'(k)*} 


(A26) 

(A27) 


It is convenient for plotting purposes and in expressing the aerodynamic coefficients to make the following definitions: 


Tt—t, — 
IT 


,Br 1 

7T 

.Br 

TT 

,v 


Up = Up' — 
7 T 


(A28) 


so that t/, o}/, t p ', and are functions of the Bm only. The variations of these four parameters with Bm are shown in 
figure 3. 

The velocity potentials for the rolling and yawing motions, completely defined by equations (A16) and (A17) and the 
constants given in equations (A21), (A24), (A25), (A26), (A27), and (A28), may now be written as 


<p v =‘px t (r v '6-\- cj p 'm) Vo(m— 0) (A29) 

and 

<p T =Br3?(j r 'd-\-u T 'm) yo(m—8) (A30) 

The pressure coefficients for the rolling and yawing motions found from equations (A29), (A30), (AlO), and (All) are 


and 


/A P\ 2pTp'm^-\-3m i (0p'xz—2mc0 r 'z l 
V 2 A V V z(mx—z ) 

(A31) 

/AP\ 2Br T/m^y-Sm^oj/xz— 2mu/z 2 
V 5 A V -jziyrvx—z) 

(A32) 
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APPENDIX B 

DEVELOPMENT OF EQUATIONS RELATING THE v - VELOCITY TO THE DISTRIBUTION FUNCTION f (<r) IN THE j/= 0 PLANE 

Equation (A8) gives tlie expression for the velocity potential (everywhere in space) resulting from a distribution of 
doublets in the xz-plane with the strength of each doublet in this distribution being governed by the doublet-strength function 
A(x,z). The derivative of this velocity potential with respect to any one of the coordinates x, y, or z will give the perturbation 
velocity in that direction. Of primary interest herein is the c-veloeity or the y-derivative of this potential 


v(x,y,z)-- 


by D {x,y,z) 

by 


(Bl) 


for points on the xz-plane. Brown and Adams in reference 15 have constructed the velocity potential in space of a distribution 
of doublets by the following approach. First, by using equations (A8) and (A18), the potential of a line of doublets in the 
xz-plane at an angle tan -I <r to the x-axis is determined. This-potential is given by 


where 


(l-BV)** 7 y 3 — 1/"^" (1-BVp 


(B2) 


X—B 2 <rZ 


7= 


V 1 — £ VVx 2 — B 2 (y s + z 2 ) 


The velocity potential of a distribution of line doublets (i. e., a surface distribution) in the xz-plane, on the vertical tail, with 
strengths governed by the distribution function /(<r) may then be written as 

<P= (B3) 

where tan -1 m is the apex angle of the vertical tail. 

Substituting equation (B2) into equation (B3) and differentiating with respect to y yield the following equation for the 
a-velocity as fiy/x approaches zero (see ref. 15) : 


v lim f r fl «-«)r By(a)Vl-W(l-B 2 < rg) i _ 3g/(<r)(l-B 2 ^)coth-y 0 2£/(<r)Vl -W 

x“e^0\Jo L (1 -£V) 2 (B<r-£0) 2 (1 -B*o 2 ) 6 '* + (1 -Wo*) 1 


C Bm r £/( <r) VI - W(1 -Rcey ZBj{ ff )0—R<je)cothr l y 0 2£/(<r)Vl 

J fl c« + .)L (1 (l-BV) 8 / 2 + (1 -w 


B 2 e 2 
£V) J J 


d{Bc) + 
2/(g)VI 


d(Bo) 


I 


(B4) 


The singularity which occurs in the - a _ - term of equation (B2) when y is set equal to zero has been accounted for in equation 
(B4). 

By taking the first and second derivatives of equation (B4) with respect to d, two other useful relations are obtained. 
They are given in the appendix of reference 15 as 


b(v/x) 

d0 


-lira 

e — *0 


{r 


*> r 3-ffV(<r) coth- 1 7o B 2 {ZBe+2Bd+B9B i o x )f(o--) _ 


BdRj(a) 


m*) 


{Bc-Bey 


L (1 -£V) 5/S 


yi-W(i-£V) 2 vi -RP&o-Bey vi - w(i ~R<R) (b*-bo) 


-f 


~ ZR<rfW) coth- 1 to B 2 (ZBr+2B6+B8B i o s )f(v) , BdBJj a) 

«(»+«> L (i-B 2 ^) 512 yi-w(i-£V ) 2 + yi ' z W{Ba-Bey 






and 


V 1 - W( 1 _BV) (Bcr—Bd) (Bv-Bdf J 


d(Bcr)- 


2B 2 frf(d) 4V1 -Wffj- 


eVl -W 




(B5) 


^= l % {evi=3¥ (B6) 

^2/ vfx ) 

The factor multiplying the f"(6) term of the expression for — as it appears in reference 15 is slightly in error and has 

been corrected in equation (B6). (The symbols f'(8) and }"{0) denote the first and second derivatives, respectively, with 
respect to 6.) 

Considering equations (A5) and (A7), it is evident that equation (B6) must be zero for both the rolling and yawing cases. 
This equation has already been satisfied by/(0) (eq. (A17)), since equation (B6) is in essence the integral equation which was 
inverted to obtain the general pressure expression from which J(d) was derived (see ref. 26). 
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APPENDIX C 

INTEGRATIONS TO OBTAIN v/x 
Tlie expression for v/x is, for 6 approaching 0 (see eq. (A19)), 


v i r r 
x~Bj o L 


— 3(TBa-\-uBm) — Bo) tank 1 -yJl—B 2 a 2 

(1 -B'a*) 6 * 


2 (tB o — |— coBm) -\JB a (Bin — Ba) 
(l-B'o*)* 


J d(Ba) 




' B " (rB<T+coBm)^Ba(Bm-B<r) (l-B^Y^l-BW Jfn , 

— (i -Bt^iBa-Bey d{B<y) 


2(rB0 J roBm) -<JB0 (Bm—Bd) Vl — 
B 2 e 




This expression has been broken into parts as indicated by the circled numbers with the third part being broken into two 
additional parts ©a and ®b because of the singularity in the integrand. Since © and ® are elementary integrations similar 
to those found in most integral tables (see ref. 27), only © will be dealt with in detail. Performing the integrations © and 
© and combining the results yield 

V1-.W f 2 (TBd+u)s/Bd(Bm-Be) ir[T(-7B6Bm-2+10Bd-Bm)+Z(4-7Bm+'iBe-BeBm)] , 

5 \ Be 8-y/l-BmO—Bd ) + 

tt[t(— 2— 10BO+Bm— 7B8Bm)-\-u>(— 4+4Z?fl— 7Bm+Bd Bm)] ^ , p9 , 

8^1+Bm(l+B0) J 

where « =wBm. The first term of expression (C2) exactly cancels ©, so that the total of ©, ©, and © for 0-*O is 

7r l~T(2-\-Bm)— co(4— iBm) , r(— 2+Bm)— w(4i-\-7Bm)~\ . 


°B |_ iJl—Bm 

The following two integrals comprising © remain to be evaluated: 


/l+Sra 


— 3 C Bm TBcr-\jBcr(Bm — I?<7)tanh. I i/l ~—B*o* 
~Bjo (1 -B 2 o*) 1 -' 2 

— 3 r Bn w-^BtriBm — Be) tank -1 Vi — BPo 2 
~B~Jo 


d(Ba ) 


(1 -B 2 ^ 


d(Ba) 


It might be mentioned at this point that the integrands of expressions (C4) and (C5) are finite and continuous over the interval 
0 to Bm and therefore must yield a finite quantity when integrated. 


The integration of expression (C4) by parts gives 


— t V-S < r (Bm Bcr ) tanh. _i yi— BV 
B(l-B*o*)* 


f Bm (Bm—2Bo) tank-yi— -BV JfT) A 
’ Jo 2 (1 -B i o*)^B<r(Bm-B<T) { * Jo 


-JBcr(Bm — Bo) 
Ba(l-B V) 2 


¥ d(Bcr) (C6) 


Integration of expression (C5) by parts gives 


— u-\IBci(,Bm—Ba ■) tanhr 1 -^/!— BPo 2 
_ jBB<r(l-SV)* _ 


' Bn wBm tank -1 -^!— B 2 * 1 £ Bm 

2BBff{l-B 2 o 2 )^^B<r(Bm^W) ® * } Jo B. 


w^Bcr(Bm — Bcr) 
BBdo* (1-BV) 2 


{ d(B<r) (07) 
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Combining expressions (C6) and (C7) results in 

^ P — (rBa+Z>) — Bo) tantr^Vl — B 3 ^ p” f (rBi7-{-it})^fWadBni—Baj 7 ,„ . . 

BB<r(l-.BV)» Jo Jo -BBV( 1 - 5 V ) 5 W)-* 


p B * [r. Bcr(Bm—2Bo-)—uBm] tanh -1 Vl— i?V 
Jo 2BBcr(l-B i a z )*-jB<T(Bm-B<r) * 


The first term of equation (C8), when evaluated at the limits, is either zero or infinity. The integrand of ®, as was noted, 
is finite over the whole interval; therefore, infini ties introduced as a result of parts integrations must, in the end, cancel 
themselves. 

The second term of equation (C8) is an elementary integration which when evaluated (with infinities neglected) yields 


n r ^2r+4co — 3 rBm — 5o)Bm — 4co — 5o)Sm-l-2r+3r5wi 

V -yJl—Bm Vl +Bm 


(C9) 


It is now convenient in integrating the third term in equation (C8) to introduce the variable substitution 


Ba 


v I'+k 

i+H 


so that Bm and k are related by 


Bm= 


2k 

1+& 2 


(CIO) 


(Oil) 


The third term in equation (C8) when transformed by equation (CIO) may be written in the form 

1 7 

± v T 

syr+F(i-p) i=i 1 

where 


(Cl 2) 


I 2 —— J* k{T+ wk?)Ftt)di 


"k( 1+A*)+35A* 
1-^ 

wWQc*- 2) 


F(^)dip 


/. 

L 


f* wk?0—k?) 

“J-»d- 
J- t a- 


WP-V) 

s*(i— *V)^ 






Fi'lddi 


tanh' 


_t ^ yi=pyi-^ 


TO)= 


l+fc* 


) 


V(i -**)(**-**) 

The integrals I 1 , I 6 , and I 7 are elementary and may be determined by an integration by parts. If the multiplicative factor 
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before the summation sign in equation (012) is neglected until all the components are totaled, these three components become 


r , r , r — rA:(l-l-A^)ir ci>A^(l+A^)ir 

i x +i 6 +i v - — p IHp 

Consider the integration required for 7 S , that is, 

df 


(C13) 


tanh 1 

V(l-A^)(l-^ 

l+AV' 

J-t v 

O—WiV-V) 


(C14) 


Let sin 0 ; then expression (C14) becomes 


X 


T/i tanh -1 • 


V(l— Ar*)(l — P sin 2 0) 
1+Fsin9 


de 


It can be shown that 


Vl — A 2 sin 2 0 

r VSEmEV^ Ltmh- timh- g £ 0S#s»/2 

L 1+Psin0 J Vl— ** sin 2 0 Vl-Jfe*sin*0 


, , Vl — A* , — Vl— ^ sin 0 

=tanh _1 ~7 === ■= +tanh~ 1 • 


Vl— t 2 sin 2 0 


Vl— i 2 sin 2 0 


0^0^— tt/2 


(015) 


(016) 


This fact allows expression (Cl5) to be written as 


tanh- (VS=r) <» „ tanh- ( M tanh- ( 

Vyi — k? sm g 0/ |' x/2 V V 1 — Ar 2 sinW _j_ J* 0 V Vl — Ar 2 sin 2 0 / 


Vl— i 2 sin 2 0 


f r/S Vl-FsinV _ r 

J-r/3 Vl— ^ sm 2 0 Jo 

The last two integrals of expression (C17) cancel each other and leave 

tanh- 1 ( ) dB 

2 f x/2 VVl-^sinV 


d0 


(Cl 7) 


«■/» Vl — Ar 2 sin 2 0 


(CIS) 


Vl — Ar* sin 2 0 

After the inverse hyperbolic tangent is replaced by its logarit hmi c equivalent and the additional variable transformation 

1-jfc 2 


sin 2 v= 


1 — Ar 2 sin 2 0 


is introduced, expression (C18) becomes 


f x/2 A 1 +sin A dv 

Jsin-tVI 3 ! 5 V ° gs 1 — sin v/ V sin 2 v— ( 1 — ¥) 


(Cl 9) 
(C20) 


It is now convenient to make the substitution 
Expression (C20) then becomes 


5=sin -i yi— A t 2 

f xt2 Xlno- 1+sin dv 

J i \ ^*1 — sin v) Vsin 2 v — si 


sin 2 5 


(C21) 


which is exactly in the form of the fourth integration formula of table 335 in reference 28. This formula gives the value of 
expression (C21) as rK'(k). The integration of 7 2 may now be expressed as 


I 2 =-k( T +uP)TrK'(k) 


(C22) 


Using the same integrating procedure for I 3 , I it and J e as just outlined for I 2 and the integration formulas in tables 335 and 
336 of reference 28 leads to 

[r*(l +AV)+35A^][1 +K' ( k ) - E' (A:)]7r 
1-Jfc* 


h 
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-5(i-4*)4**Kr(jfc)- 


5*M1 +K'(k)-E'(k)] 
l-k* 


where n (1 — A 2 ) , Vl — A*> U is a complete elliptic integral of the third kind with modulus -Jl—ic* and parameter — (1 —Ic 2 ). 

S ummin g all the various parts contributing to the third term in equation (C8), including the common factor, gives the 
following expression: 


-1 f , xtf / (*)[2**T+SA*(l+**)] 

5o=^CTP\ +1 + r=F 

* E ' (k) [2k^+rk(l + P )j + - Fx(1 _ P)n (i i |J J- (C23) 


The addition of expression (C23) to expression (C9) completely evaluates Q. Expression (C3) gives the evaluation of ®, ©, 
and ©. Before writing the total integration, that is the sum of expressions (C23), (C9), and (C3), it is desirable to combine 
•expressions (C3) and (C9), which are functions of Bm, and transform them by equation (Cll) to functions of k. This pro- 


cedure yields 


~ (ZiAr 2 — (— rA 3 ) 
B^jT+V{l-k*) 


(C24) 


The total integration may now be written in terms of the parameter k as 

^==^-^|iC'(A)[2Fr+SF(l+A I )]-S / (A)[2A i 5+rA(l+A 2 )]+5A s (l-A 3 ) 3 n|^-(l-A 2 ),l-A 2 ,|Jj. (025) 

By use of the process commonly known as interchanging the amplitude and parameter (see pp. 133 to 141 of ref. 29) the elliptic 
integral of the third kind appearing in equation (C25) is found to be equivalent to — p — This operation permits the expres- 
sion for v/x (eq. (C25)) to assume the form given in equation (A20). 
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TABLE I 

FORMULAS FOR POTENTIAL DISTRIBUTION DUE TO CONSTANT SIDESLIP 



Region (see 
sketch) 

?(x,0 + ,r) 

I 

z) 

-jBhn'-X 

II 

V P r frrr -) co;r > mx-z(2Bm-D 

injBhrt-lL nx-* 

2 V zm(x— Bz) ( Bm — 1) J 

III 

V P f T ) co _, «^-z+2(l + Ba»)(z-6) 

irjBBn'- mx ~ z 

2 y(z—b) (Bm+1)[ XI +#*»)— m{x-\-Bz)] > 

' 

IV | 2 (11+ III -I) 

V 

Vf}H(Bm)^z{mx — z) 
Bm 
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TABLE IK 

FORMULAS FOR POTENTIAL DISTRIBUTION DUE TO STEADY ROLLING 



,>Modi lines 



(Tail is in xs-plane; positive rolling) 


Region 
(see sketch) 




V 

2 (j 8 *m» 


{ ( mx ~ 2 )[ — mx+z(2B*m*— 1)]| 


>-!)*/»{ 


mi(4Bm— l) + s(2g a nt > — 2Bm— 3) lmz(x—Bz) 


3 (Bhn 3 — 1) V Bm+l 

<jnx—z)[mx—z(2B t m l —l)\ mx+e(l — 2Bm) ) 

2(B*»»«-l)»/» 008 } 


(mx— z)[— mx+z(2B*m s — 1)] mi— z+2{l+Bm)(z— b) 

2 cos 

mx—z 


(mi— 6) (4Bm— 1) + (z— 6) (3 — 2B t nP—2Bm) — 6b(B 2 m} — 1) „ — — — ) 

^ ^ - ^V(2-6)(B>n+l)[6(l+Bm)-m(i+Br)]J 


z(ii+in-i) 


V(r v / z-\-<a v , mx)^z(mx—z) 
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TABLE IV 

FORMULAS FOR PRESSURE-DIFFERENCE COEFFICIENT DUE TO STEADY ROLLING 



Region 
(see sketch) 

&p . . 

— (x,z) 

' E9 1 

— mx) 

| H 

V(BW- 1 )** 

ii 

4 p F 2 Bm 1 jzm(x—Bz) 

xV LCR^’-l) V Bm+1 


m t {x—B t zm) , mx+z(l — 2Bm)~| 

(Bhn'—V)*! 1 COS_ mx-z J 

hi 

4 pm , . , tra-i+2(l+Bm)(z-li) 

rr/n. .. — rrrs 1 ( Bhn'z—mx ) cost 1 — — 

1) 3/5 ( wrc— z 

) 


2Bm V (z — 6) (Bm + 1 ) (6 ( 1 + Bm) — m {x + Bz) ] | 

IV 

S (11 + III -I) 


2 p r p ' mz ! + 3m 5 u J ,'zz — 2»iMp' z 1 


V 


Vz(mx— z) 
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TABLE V 

FORMULAS FOR POTENTIAL DISTRIBUTION DUE TO STEADY YAWING 



Region 
(ace sketch) 

< p ( x , 0 + , z ) 

I 

[ ^ ( 2m + B * m *) + 2xz B'z'm] 

II 

r (x(5B t m i +4Bm—Q)—B L mz(2Bm+l) Imz(x—Bz) 

x { 3 (13*771* — 1) V Bto+1 1 

(mx—z)[mzB 1 +x(m t B 2 — 2)] , mx— z(2Bm— 1) ) 

2(J3»m*-l)*'» C0S mx-z ] 

m 

r$mx(5Bhn*-‘lBm-fi)+B*m'(z-b)(2Bm-l)+Bmb(4:+Bm) r tttttt ■ B ) , . n j . 

(mx-z)[B 1 m*(nix-\-z)—2mx] , mx— z+2{l-\-Bm) (z—b) ) 

2m(Bhv?—\y li C0S ~ mx-z ) 

IV 

2(11 + 111-1) 

V 

Br (t,' z + w/ mx) V z (.mx — z) 
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TABLE VI 

FORMULAS FOR PRESSURE-DIFFERENCE COEFFICIENT DUE TO STEADY YAWING 
U / I 



(Tail is in xz-plane; positive yawing) 


Region 
(see sketch) 

aP , . 

— (*,*) 

I 

V{Bhn*-\yn 2 ) + «] 

II 

;vw*=w* 

. ) 


2(B 5 ot s -)-Bto— 1) V mz(Bm — l)(x— Bz) > 

m 

4r f r. , om , mx- z+2(l+Bm)(e-b) , 

xV(B»m 1 -l)*fl + 2)] cos" + 


2(B 2 m t —Bm—l) ->J(z—b)(Bm+l)[b(l+Bm)—m(x+Bz)] > 

IV 

sqi+in-i) 

V 

2Br r r ' mz 5 + 3m 2 w/:rz — 2mw/ z 1 


z(mx—z) 
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Region 
(see sketch) 


9>(x,0+, z) 


P ( mx—z ) 


-JBhn'-l 


-ivsSrdv^-* >(B ”-D [w + «^=|g i |^ a ±2 ] + 

[— ”~f fr 11 ] [iSSgff +«-] j 


5={V(*-6)(B*+l)[6(l+5m)-w(*+B*)] [2tF+ AP(4+jB ” t)(6 & ) C 2 g»i— l) ~j + 


m- z+2(1+Bot)(z— 6)' 


Mhniz—mz) 
' 1) 


•J {mx—z) | 


s(ii+m-D 
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TABLE VIII 

FORMULAS FOR PRESSURE-DIFFERENCE COEFFICIENT DUE TO CONSTANT LATERAL ACCELERATION 



(Tail is in 12-plane; positive fi in positive ^-direction) 


Region 
(see sketch) 

AP, . 

I 

4 f}(mx—z) (l + n*) 
V(B*m 5 - 1 )M* 

II 

40 |B *(1 + m’)(rax— z) , mi-i( 2 Bm-l) 

it&Vt/BW- 1 t B , m t — 1 mx-z 

2 [l-— ( ^t 5 f—^]v^(i-gz)(R ro -l)} 

III 

4/3 |B s (l + m I )(77i2:— z) mx—z + 2 (l+Bm)(z—b) 

r&V-J&ri 1 — 1 t Bhn 1 — 1 C03 mx—z 

MRB^-Bm-in , 1 

L BVn’ — 1 J V(*~ &KBi»+l)[6(l+B»t) — jn(x+Bz)]J- 

IV 

2(11 + 111-1) 

V 

20 r(B=+ 1)B(t/77iz j + 3 otV'iz— 2W* 1 ) 

JS’H ->Jz{mx—z) 

(B 1 +l)xzB'(Bm) 2 H(JBm)-<] z(mx— z)”| 

B-<Jz(rnx—z) Bm J 
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TABLE IX 

STABILITY DERIVATIVES FOR HALF-DELTA AND RECTANGULAR ISOLATED VERTIC AL TAILS 


Derivative 

(a) 

Half-delta tails 

Rectangular tails 
ABil 

0^AB<2 

ABg2 


~ H(Bm) 

4 / AB 

Fvab+2 

B L 2AB/ 

Cn n 


16 / AB 

3ABVAB+2 

— P-- M 

AB V 3AB/ 

S 

-Wi H ^ 

4 AB+1 

3S VAB(AB+ 2) 

B V 2AB J 


(V + 2^0 

4 VAB(AB+3) 
3B (AB+2)" 2 

2AB-1 

AB 2 


-y (V + 20,/) 

2(AB+3) 

VAB(AB+2)»« 

3AB-2 

3AW 


tAB / 5 , , A 

45 \8 T ” +<Jp ) 

2A 2 B 2 +6AB+3 

3BVAB(AB+2)w» 

1+ 4AB - 24A 2 B 2 + 32 A 3 B» 
24A 1 B< 

Cr r 

^ /■> 

8(2AB+6) 

2(3AB — 1) 


3VAB(AB+2)w 

3 A 2 B> 

C'r 

_3*B . / , 

4B(2AB+5) 

B(8AB— 3) 


IAB(AB+2)P« 

6A»B S 

C‘ r 

*AB ( 5 A 

4 U Tr+ “ r J 

2(3A 2 B>+9AB+5) 

3[AB(AB+2)P rt 

3AB-1 
3 A*B» 

Cr i 

ABt r , 

4 L 3 A’B 2 J 

ABt r , 8 H(Bm)-l 

4 B 2 L + A 2 B* J 

8 B 1 — AB— 1 

353 VAB(AB+ 2)* 

2(B 1 +2-r3AB) 

3A’B< 

s 

4 L2 Tr+3 “ r A’B 2 J + 

t T3 , , „ , 12/f(Bm)"l 
4B La r + 3 r A*B* J 

4 B 2 — AB— 1 

5 VA»B3(AB+2)’ 

3B 2 — 8AB+ 6 
6A*B* 


8 L 2 ‘ 4 Tr A 2 B' J 

ABt T 5 , 8H(Bm)~\ 

8 B> L r 4 r A 2 B« J 

2 B 1 +A 2 B>+3AB+3 
3S 5 VA , B 2 (AB+2)» 

B 2 +2-3AB 

3A J B* 


* Angular velocities and moments measured about the system of body axes shown in figure 2 (b). 


435870 — 57 27 
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TABLE X 

STABILITY DERIVATIVES FOR TRIANGULAR ISOLATED VERTICAL TAILS 
WITH SUBSONIC LEADING EDGE AND SUPERSONIC TRAILING EDGE 


Derivative 

(a) 

Formula 

(b) 


-iVa 

G* 

x(3AB+2Bm) 

GBmyjZBmAB 

C>, 


C r, 

fV2 BmAB^+^+^w) 

Cn p 

■nj2BmAB £t p ( 32 Bm + 16AB) + “ 1 ' (.S2Bm + 4AB 1 8A»bO] 

Ci, 

t ^V2B«ab(|v+| 

Cr r 

W2BmAB(^+^+^ov') 

Cn T 

xBV2 BmAB [r r ( 32Bm +i6AB) + “ r (32 Bm h 4 AB 1 8A*B>)] 

Cl r 

^WwAbQt/+^'+^ «/) 

Cr i 

,„[>/ , «//, , 2Bm\ II (Bm) H(Bm)l 
in/2BmAB l^ + ABj AB*m* GBmABj 

«■ )-a- 1 f? r Tr ' 1 ^ Y 1 1 *<£») (-Bm) 1 

^ V2BmA5 L4 + 4 (,1+AB; 4B*m’ 2BmABj 

s 

xBV2BmAB^r r ( 32 Bm + 16AB) + " r (s2Bm ^4AB 1 8 A’B 1 ) 

(32B3m l+ 8B 5 m J AB + 8BmA s B 1 ) + 

B V2BmAB [\> ( 32Bjn 1 16AB ) 1 (32Bm + 4AB + 8A*B J ) 

(32B’m» + 4B 3 m 2 AB + 8BmA J B ! ) 


3 T 2 V2BmAB[5r r '+ 5u / + 6 f-./ (J m ,+ EmAs) 
g^VMB^'+fiov'+^V (r^+R A*) *(*">] 


* Angular velocities and moments measured about the system of body axes shown in figure 2 (b) . 

b Formulas valid for triangular plan forms with either sweptback or sweptforward trailing edge provided the edge is supersonic. 
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TABLE XI 

TRAN8FER-OF-AXES FORMULAS 

Stability derivatives Formulas for transfer to a body system 
in a body system of of axes with origin displaced dis- 

axes with origin at tances xa (positive forward) and zo 
tail apex (see fig. 2 (positive downward) from the tail 

(b)) apex (see fig. 2 (c)) 
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\c. 



(b) 


' (a) Plan forms of vertical tails analyzed. (Trailing edge may be 
either sweptbaok or sweptforward provided it remains supersonic.) 

(b) Special cases (rectangular and half-delta vertical tails) for which 
stability-derivative curves are presented. 

Figure 1 . — Tail plan forms and associated data. 



(a) Body-axes system used for analysis. Free-stream veloolty V. 

(b) Principal body-axes system used for presentation of stability 
derivatives. Entire system moving with flight velocity V. 

(c) Same type of axes system as (b) with origin translated. 

Figure 2. — Systems of body axes. Positive directions of axes, forces, 
and moments indicated by arrows. 




Bmur, Qr| d 



-BC,., or C„ 
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Fioube 5 . — Curves for determining the stability derivatives due to steady rolling Cr v , C„ p , and C, p f or isolated half-delta vertical tails. Deriva- 
tives baaed on vertical-tail parameters b, S, and angle pb/ V; principal body-axes system -with origin at leading edge of root section. 
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Figure 6. — Curves for determining the stability derivatives due to steady yawing Ct t , C, r , and Ci r for isolated half-delta vertical tails. Deriva- 
tives based on vertical-tail parameters b, S, and angle rbjV; principal body-axes system with origin at leading edge of root seotlon. 
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FiatjRB 7. — Curves for determining the stability derivative due to constant lateral acceleration Or- for isolated half-delta vertical tails. Deriva- 

fi 

tive based on vertical-tali parameters b, S, and angle $b/V; principal body-axes system. 




430875 — 57 28 
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Figube 8. — Curves for determining the stability derivative due to 
tive based on vertical-tail parameters ft, S, and angle $b/V; 


constant lateral acceleration C n . for isolated half-delta vertioal tails. Derive- 
principal body-axes system with origin at leading edge of root section. 
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Fioueb 9. — Curves for determining the stability derivative due to constant lateral acceleration Cj. for isolated half-delta vertical tails. Derive, 
tive based on vertical-tail parameters 6, S, and angle $b/V; principal body-axes system with origin at leading edge of root section. 
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Fxqueb 10. — Curves for determining the stability derivatives due to constant sideslip Cr P C n p and Ci ff for isolated reotangular vertical tails. 
Derivatives based on vertical-tail parameters b and S; principal body-axes system with origin at leading edge of root seotion. 
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Figukb 11. — Curves for determining the stability derivatives due to steady rolling CV^, C Bj>) and Ci v for isolated rectangular vertical tails. 
Derivatives based on vertical-tail parameters b, S, and angle pb/V; principal body-axes system with origin at leading edge of root section. 




420 


REPORT 1268 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Floras 12. — Curves for determining the stability derivatives due to steady yawing Cr T , C, r , and Ci r for isolated rectangular vertical tails. Deriva- 
tives based on vertical-tail parameters 6, S, and angle rb/V; principal body-axes system with origin at leading edge of root section. 
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Fiqdbb 13. — Curves for determining the stability derivatives due to constant lateral acceleration Cy- and Cj. for isolated rectangular vertical 

& s 

tails. Derivatives based on vertical-tail parameters b, S, and angle $b/V; principal body-axes system with origin at leading edgo of root 
section. 


i 
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Fiqube 14. — Curves for determining the stability derivative due to constant lateral acceleration C.. for isolated rectangular vertical tails. 

fi 

Derivative based on vertical-tail parameters b, S, and angle $b/V; principal body-axes system with origin at leading edge of root seotion. 
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(b) Steady yawing. 
Figttbe 15. — Continued. 









(c) Constant lateral acceleration. 
Figubb 15. — Concluded. 
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